Statefinder hierarchy of bimetric and galileon models for concordance cosmology 
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In this paper, we use statefinder hierarchy method to distinguish between bimetric theory of 
massive gravity and DGP model applied to late time expansion of universe. We also carry out com- 
parison of bimetric model and galileon system . We show that statefinder diagnostic can differentiate 
between the currently discussed cosmological models. 
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I. INTRODUCTION 



There exists a belief that observed late time cosmic 
acceleration is driven by some unknown exotic energy 
component characterized by negative pressure dubbed 
1 dark energy' [3-|3|- This hypothesis is supported by 
number of observational results related to Type la su- 
pernovae [l], Q , cosmic microwave background radiation 
and large scale structures 

The simplest candidate of dark energy is cosmological 
constant A with p = —p. There is also a variety 
of dynamical dark energy models which can fit the 
observations. In view of the forthcoming observations, 
it is of utmost importance to find ways to distinguish 
these models. Different diagnostic measures have been 
proposed in the literature to distinguish dark energy 
models; Om diagnostic @ and statefinder [f| 0] are two 
examples for such diagnostic measures (see also Ref. [IH 
on the related theme). In this paper we shall employ 
statefinder for distinguishing some recently proposed 
models of dark energy The statefinder pair {r, s} is a 
geometrical diagnostic which is algebraically related to 
the higher derivative of scale factor "a" with respect 
to time. The deceleration parameter (q), statefinder 
(r) and snap (0:4) contains second, third and fourth 
order derivative of scale factor respectively. It is really 
interesting that statefinder can successfully differentiate 
between a large variety of models 0, Q ■ 

The statefinder diagnostics for nonminimally coupled 
scalar field system and galileon field which is generically 
nonminimally coupled system has been investigated in 
Ref. In which we have shown that the nonmini- 

mally coupled scalar field and galileon models are suc- 
cessfully differentiated from other popular dark energy 
models such as chaplygin gas, quintessence and Dvali, 
Gabadadze and Porrati (DGP) models in r - s, r - q plane. 
Recently, a more refined diagnostic known as 1 Statefinder 
hierarchy is introduced in Ref. ■ 

Recently, statefinder pair {r, s} was used to compare 
nonminimally coupled scalar field and galileon models 
with other popular dark energy models We have 



shown that statefinder pairs successfully differentiate 
between LCDM, nonminimal, galileon, chaplygin gas, 
quintessence, DGP on the r - s, r - q plane. 

In the present paper we consider bimetric model of 
massive gravity which is closely related to nonlinear mas- 
sive gravity a la dRGT [l^ - [f^ | . An interesting scalar field 
dubbed galileon despite the higher order derivative terms 
does not suffer from Ostrogradki's ghosts. Galilcons 
emerge in dRGT in the so called decoupling limit which 
is a valid limit for studying the Vainshtein mechanism 
and galileon is a natural device to implement the latter. 
In fact, the lower order galileon lagrangian is responsible 
for the consistency of DGP with local physics. A large 
number of papers are devoted to cosmological dynamics 
of galileon field. It was first demonstrated in Ref. [f5[ 
that one needs a higher order galileon system to achieve 
de Sitter solution. 

It is interesting to distinguish the massive gravity mod- 
els from models based on galileon field. In this paper, we 
use statefinder pair {r, s} to differentiate between bimet- 
ric and galileon models; we also use statefinder hierarchy 
S n to compare bimetric model with other popular dark 
energy models. 



II. THE STATEFINDER HIERARCHY 

In what follows we shall work in the framework of spa- 
tially homogeneous and isotropic universe; in this case, 
the scale factor a(t) is the only dynamical variable. Since 
we shall be interested in the late time behavior of expan- 
sion of universe, we consider the Taylor expansion of the 
scale factor around the present epoch to Gil : 

a(t) = o(<o) + a(t Q ) V W - *»)]" ■ (!) 
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The deceleration parameter is defined as 
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The statefinder pair {r, s} and the Snap are defined as 
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The Statefinder hierarchy 5„ is given by [Toj : 
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and so on, 
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where 012, 03, a4 are given by ((3]), |4]), (j6|) and 

ft 0m (l+z) 3 /^). 

The alternate form of the Statefinder 
S3 := S 3 , 

Si 1] := ai + 3{l + q), 
The second Statefinder corresponding to (j9|), 
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therefore, {S^^S^} = {1,0}. In the discussion to fol- 
low, we shall use the above defined parameters to distin- 
guish various models of dark energy. 



III. DARK ENERGY MODELS 

• Bimetric theory of massive gravity following 
Refs. [H E3 we consider the bimetric massive 
gravity action 
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S = / d 4 x^dctgR(g) - -J- / d 4 x^detfR(f) 
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FIG. 1: Evolution of u> e // versus redshift z for bimetric 
model of massive gravity. 

and / metric is a reference metric. This theory is 
ghost free and reproduces dRGT in a certain limit. 

The generalized Friedmann equations for flat Uni- 
verse in the bimetric theory are given by (k = 0) 
10: 
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This theory contains two space-time metrics g and 
/. The g metric is assumed to be the physical metric 



0, 
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where, p m and p r are the energy density of matter 
and radiation respectively. 

The Hubble parameter and w e ff for this model 
have the form 
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FIG. 2: The left panel (a) shows the time evolution of the statefinder pair {r, s} for bimetric (silod line) and DGP (dashed 
line) models. Bimetric model lies to the left of the LCDM fixed point (r = 1, s = 0) and DGP model lies to the right of the 
LCDM fixed point (r = 1, s = 0). For Bimetric model, s remains at zero, whereas r first increases from unity to a maximum 
value, then decreases to unity. For DGP model, s decreases to zero, whereas r first decreases from unity to a minimum value, 
then increases to unity. Both models converge to the fixed point (r = 1, s = 0) which corresponds to LCDM. The right panel 
(b) shows the time evolution of the statefinder pair {r, q} for Bimetric (silod line) and DGP (dashed line) models. Both models 
diverge at the same point (r = l,q = 0.5) which corresponds to a matter dominated universe (SCDM), and converge to the 
same point (r = 1, q = —1) which corresponds to the de Sitter expansion(dS). The dark dots on the curves show current values 
{?"o,so} (left) and {ro,go} (right). In both models, f2o m = 0.3 at the current epoch. 
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where, Bi = m 2 [3i/Ho, m? is the order of H$ and 
Pi is absorbed into m 2 ; y = Y/a, Y and a are the 
scale factors corresponding to the metric / and g 
respectively. Here ' denotes derivative with respect 
to Ina. 



• Galileon model (15H1 

The Galileon is a massless scalar field (it); whose 
action is invariant under the Galilean transforma- 
tion ir(x) — > 7t(.t) + b^x^ + c, where and c are 
a constant four vector and scalar rspectively. The 
galileon action is of the form [l5| 
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term which occurs in DGP; L4 and L5 are higher 
order Lagrangians. In this case, the evolution equa- 
tions in a spatially flat background have the form 
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- —IT - C 3 TT 7T 

-^c 4 tt 3 (3H 2 ir + 2Htt + 8H^j , (17) 



(ip m = -c 2 (3Htt + + 3c 3 n (3H 2 ir + Hn + 2Hj^J 
-18c 4 i?7T 2 (iH 2 i: + 2Htt + 3Hn) . (18) 



where, {ct} are constants, /3 is coupling of field with 
matter and L" are Lagrangians for galileon 
field. Li is linear in field and is often omitted as- 
suming ci =0, L 2 represents the standard kinetic 
term, L 3 = ((9 M 7r) 2 D7r is the third order galileon 



The Hubble parameter for galileon model has the 
form 

H(z) = H b 0m (l + zf + + zf^+^A V2 , (19) 
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FIG. 3: This figure shows the evolution of the Statefinder Si, S 1 ^ 1 ' plotted against fi m = f2om(l + z) 3 /h 2 . Large values 
fim — > 1 correspond to the distant past (z > 1), while small values Q m — > correspond to the future (z — > —1). The models 
are: bimetric (red), galileon (magenta), DGP (blue). The horizontal black line shows LCDM. The vertical band centered at 
firjm = 0.3 roughly corresponds to the present epoch. 
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FIG. 4: This figure shows the evolution of the Statefinder 
plotted against fl m . 
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• The DGP model hj: 
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In order to apply the statefinder analysis to the mod- 
els under consideration, we notice that the deceleration 
parameter q, statefinder pair {r, s} and Snap can be 
easily expressed in terms of Hubble parameter H(z) and 
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FIG. 5: The panels (a), (b) show the Statefinder Sg 2 \ plotted against S7 m .The dark energy models are: bimetric (red), 
galileon (magenta), DGP (blue). The horizontal black line shows LCDM. The vertical band centered at f2om = 0.3 roughly 
corresponds to the present epoch. 



its derivatives as follows: 
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In Fig. 21 galileon model is highly degenerate, DGP 
model is nearly degenerate whereas bimetric model of 
massive gravity is showing nondcgencracy around the 

(2) 

present value of S2 m . In S3 — TO plane, galileon model 
is nearly degenerate whereas bimetric and DGP models 
arc nondcgcncratc around the present value of f2 m see 
FigJS] (a). In S\ ' — Q m plane, galileon model is highly 
degenerate, DGP model is nearly degenerate whereas bi- 
metric model of massive gravity is nondegenerate around 
the present value of Q m see Fig[5] (b). In Fig. HI we 



exhibit the phase portrait in — S^ 1 ' plane where de- 
generacies among bimetric, galileon and DGP models are 
broken. 



,(D 



where, H is given by (flijl. (fT9|). (|20|) for different dark 
energy models. 

We are now in position to present the numerical results 
for statefinder parameters for each model and strike a 
comparision between them. In Fig. [5J We show the 
time evolution of the statefinder pairs {r, s} and {r, q} 
for bimetric model of massive gravity and DGP model. 
One can see that both the models are differentiated by 
statefinder pairs. The statefinder hierarchy is used for 
higher derivative of scale factor; which can easily distin- 
guish between LCDM, bimetric, galileon and other dark 
energy models. In Fig. |3l we show the behavior of bimet- 
ric, galileon and DGP models in S% — 0, m and S3 — fi m 
plane. All models considered in this paper show nonde- 
generacy around the present value of Q m = 0.3 in S2 — f2 m 
and So 1 ^ — fl m plane. 



Particularly, we should comment on S2 — f2 m and 

S^ 1} - fl m plane. From Fig. H one can clearly see that 

S2 and S3 1 ' does not have degeneracies between bimet- 
ric, galileon and DGP models around the present value of 



S3 plane degeneracies arc bro- 



O m . Similarly, in S 4 
ken among bimetric, galileon and DGP models around 
the present value of tt m see Fig. [SI In all figures we con- 
sider f2om = 0.3, for other values of f2o m the conclusions 
also remain same. Other statefinder hierarchy does not 
seem to perform well in distinguishing the different dark 
energy models considered in this paper. 
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FIG. 6: This figure shows the Statefinders and S3 
for bimetric (red), galileon (magenta), DGP (blue) models. 
LCDM corresponds to fixed point (1,1). The arrows and dots 
show time evolution and present epoch in different models 
respectively, we assumed fiom = 0.3 (present epoch). 



by using the statefinder pairs {r, s} and {r, q}. We also 
carried out comparison between bimetric theory of mas- 
sive gravity, galileon modified gravity and other popu- 
lar dark energy models using the statefinder hierarchy in 
concordance cosmology Our investigation in S2 — O m 
and — fl m plane shows no degeneracy among bimet- 
ric, galileon and DGP models around the present value of 
f2 m and all models considered in this paper are success- 
fully differentiated by statefinder hierarchy on this plane. 
We have also noticed that our analysis presents the best 
comparison among the different dark energy models con- 
sidered in this paper. Figures ([3l [6]) show nondegeneracy 
between bimetric, galileon and other popular dark en- 
ergy models around the present value of fl m . We find 
that S^', and does not perform well in distin- 
guishing the different dark energy models considered in 
this paper see Fig. |4] and [5] . Looking at the success of 
the statefinder hierarchy diagnostic, we are encouraged to 
apply the analysis to models of extended massive gravity 
discussed in the literature. 
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